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1. Introduction 

Mass loss phenomena in the form of powerful bipolar jets 
and molecular outflows are often associated with the early 
evolution of protostars. Together with indirectly regulat- 
ing the accretion process, they play a crucial role also in 
the interaction between the protostar and the natal envi- 
ronment, causing injection of momentum, kinetic energy 
and turbulence in the ISM, along with irreversible modi- 
fications of its chemical structure and physical conditions. 
The most violent interaction occurs at the terminal working 
surface, where the supersonic flow impacts the undisturbed 
medium and most of the ambient material is entrained by 
the jet (e.g. Reipurth & Bally, 2001). In a schematic model, 
the interaction occurs via two shocks: an internal working 
surface (Mach disk) which decelerates the jet gas and a for- 
ward shock which accelerates the ambient gas producing a 
mixture of shock velocities (e.g. Hartigan, 1989). This lat- 
ter has often a curve-shaped morphology: therefore, only 
the component of motion along the jet axis is slowed going 
from the apex of the shock toward the receding parts, where 
large transverse motions occur. In a classical scenario, for 
sufficiently fast shocks, the head of the bow is a pure dis- 
sociative (J)ump-type shock (Hollenbach & McKee, 1989), 
which changes along the bow flanks, where the impact oc- 
curs obliquely, to slower, non-dissociative (C)ontinuos-type 
shocks (Draine, 1980). In this framework, highly excitated 
ionic emission should arise at the head of the bow, while 
molecular emission mainly originates from the cooling re- 
gions along the flanks and behind the bow. Recent mod- 
els, however, predict that mixtures of J- and C-type shocks 
can occur along the overall structure of the bow, or, in the 
presence of a sufficiently strong magnetic field, that a J- 
type shock can evolve into a C-type, remaining embedded 
at early time (Smith & Rosen, 2003; Smith & Mac Low, 
1997; Flower et al., 2003; Le Bourlot et al., 2002, Lesaffre 
et al., 2004a,b). 

The large number of parameters which regulate the 
shock physics (e.g. strength and direction of the local 
magnetic field, ionisation fraction and density gradient 
between the shocked gas and the local environment) 
can be constrained only through dedicated observations 
aimed at probing the physical and kinematical properties 
of the gas all along the bow structure. In this respect, 
observations in the near-infrared represent a well suited 
tool: both H2 ro- vibrational lines, which are the main 
gas coolants in continuos shock components, and fine 
structure transitions of abundant atomic species (e.g. 
Fell, Ci, Sn), which are expected to emit in dissociative 
shocks, fall in this wavelength range. Indeed, in recent 
years multiple observations of jets and bow shocks have 
been conducted in the near-infrared by means of long-slit 
spectroscopy (e.g. Eisloffel, Smith & Davis, 2000; Giannini 
et al., 2004; Smith, Froebrich & Eisloffel, 2003; Nisini et 
al., 2002; O'Connell, Smith & Davis, 2004; Nisini et al., 
2005). The observed spectra are typically rich in both 
molecular (H2) and ionic ([Fell]) lines, implying that 
strong gradients in the local conditions (e.g. temperature, 
fractional ionisation) do occur and that multiple shock 
components are simultaneously at work. 
The main limitation of long-slit spectroscopy for studying 
in detail bow-shock morphologies arises from the poor 
coverage of the extended bow surface usually obtain- 
able in a reasonable amount of observing time; on the 



contrary, Integral Field Spectroscopy (IFU) represents 
a well-tailored tool to overcome this problem, since it 
allows us to obtain simultaneously 2-D maps at different 
wavelengths. In this paper we present the spectral images 
obtained with the IFU facility SINFONI (Eisenhauer et 
al., 2003, Bonnet et al., 2004) of a prototype bow-shock, 
namely the Herbig-Haro object HH99B. This is located 
in the RCrA molecular core at d^ 130 pc (Marraco & 
Rydgren, 1981), and it was firstly discovered in the optical 
by Hartigan & Graham (1987), who suggested it is the 
red-shifted lobe of the outflow powered by the HH100-IR 
source. More recently, Wilking et al.(1997) have proposed, 
on the basis of near-infrared images, the infrared source 
IRS9 and the Herbig Ae star RCrA as other possible 
exciting source candidates. HH99B was firstly imaged in 
the near-infrared by Davis et al. (1999, hereafter D99) 
who identified three different emission zones: one at the 
head of the bow (subsequently named B0 by M c Coey et 
al., 2004, hereafter MC04) where the bulk of the emission 
comes from ionised gas, and two bow-flanks (knots Bl 
and B3), which emit mainly in H2 lines. A further H2 
knot, immediately behind the bow apex, was identified 
as knot B2. In the framework of both bow- (D99) and 
planar- (MC04) shock models, some attempt has been 
made to model the line emission observed in HH99B: the 
H2 morphology is well fitted by a C-type bow shock, while 
the physical conditions of the molecular gas (measured 
by means of long-slit spectroscopy) have been accounted 
for by a planar J-type shock with a magnetic precursor. 
None of these two models, however, is able to reproduce 
the copious ionic emission, that requires the presence of a 
further fully dissociative shock component. 
With the present work we aim at putting strong observa- 
tional constraints on bow-shock models, (whose detailed 
application will be subject of a separated paper). In 
particular we intend to: (i) morphologically characterise 
the emission of the different lines; (i) derive maps of the 
main physical parameters which govern the shock physics; 
(in) study the velocity field(s) along the bow structure. 

Our work is organised as follows: Section 2 describes our 
observations and the obtained results; the line excitation 
analysis and the kinematical properties are then presented 
in Sections 3 and 4. Concluding remarks are given in Section 
5, while in Appendix A we describe the procedure we have 
applied to derive the inclination angle of HH99B. 



2. Observations and Results 

HH99B was observed during four different runs in May 
and July 2006 with the SINFONI-SPIFFI instrument 
at the VLT-UT4 (ESO Paranal, Chile). The coordi- 
nates of the pointed position are: a2ooo = 19 ,l 02 m 05.4 ;s , 
^20oo=-36°54'39". The Integral Field Unit was employed 
to obtain spectroscopic data in J (1.10-1.40 /mi), H 
(1.45-1.85 /im) and K (1.95-2.45 /mi) bands, at spec- 
tral resolution 2000, 3000, and 4000, respectively, each 
wavelength band fitting on the 2048 pixels of the Hawaii 
detector in the dispersion direction. An image sheer con- 
verts the bidimcnsional ficld-of-view into a one-dimensional 
long-slit, that is fed into a spectrograph to disperse the 
light of each pixel simultaneously. As a consequence seeing 
and atmospheric response can affect neither the emission 
morphology at a given wavelength nor intensity ratios of 
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lines within the same filter. The spatial resolution was 
selected at 0.25" per image slice, which corresponds to a 
field-of-view of 8"x8". No adaptive optics is supported in 
this configuration. The total integration time is 2400 s, 
f500 s, and 1800 s in J, H, and K bands, respectively. 
The observations were acquired by nodding the telescope 
in the usual ABB'A' mode and a telluric B-type standard 
star close to the source was observed to remove the 
atmospheric spectral response. 

The S1NFONI data reduction pipeline (Modigliani, 
Ballester & Peron, 2007) has been used to subtract the 
sky emission, to construct dark and bad pixel maps and 
flat field images, to correct for optical distorsions and to 
measure wavelength calibration by means of Xenon- Argon 
lamp images. Further analysis has been carried out with 
IRAF packages and IDL scripts, that were used to remove 
telluric absorption features: this task was accomplished by 
dividing the target images by those of the telluric stan- 
dard star, once corrected for both the stellar continuum 
shape (a black-body function at the stellar temperature) 
and its intrinsic absorption features (mainly hydrogen re- 
combination lines). Photometric calibration was obtained 
from the same standard star. Since there is not overlap- 
ping in the spectral range covered by the different grisms, 
no cross-calibration has been performed at this step of the 
data reduction. However, we have checked a posteriori the 
reliability of our absolute fluxes once an extinction map for 
the H2 line emission has been obtained (Sect. 3. 1.2): we have 
measured the departure from the theoretical value in the 
image ratio of de- reddened lines coming from the same level 
and lying in different spectral bands. In this way we esti- 
mate a cross-calibration better than 5% between J and H 
bands, and than 12% between H and K bands. Atmospheric 
OH lines were also used to refine wavelength calibration 
and to measure the effective spectral resolution: we obtain 
R ~ 1900, 2700, 3500 in J, H, and K bands, that corre- 
spond to 160, 110, 85 km s^ 1 , respectively. At these spec- 
tral resolutions we are able to resolve the brightest lines 
(e.g. H2 2.122/^m and [Fen] 1.644/^m), which are observed 
at S/N ratios of ~ 10 2 -10 3 . As a final step, images of lines 
observed in individual bands ([Fen]1.257/im and 1.644/im 
in the J and H bands and H2 2.122/im in the K band) were 
used to re-align the images acquired in different observing 
runs. 

As a result, we obtained a 3D data-cube containing the 
HH99B image in more than 170 linesQ As an example of 
our results, we have integrated the signal in the areas cor- 
responding to knots B0 and B3 (see FigureQJ: the corre- 
sponding spectra are shown in Figures[5]H] The large ma- 
jority of the detected lines are H2 ro-vibrational lines (121). 
For these, we list in Table [1] spectral identification, vac- 
uum wavelength, excitation energy (in K) and the maxi- 
mum S/N ratio registered in the corresponding image. The 
detected ro-vibrational transitions come from levels with v 
< 7 and E up up to ~ 38 000 K, many of them never ob- 
served before in HH objects. In particular, as we show in 
Figure [1] (upper panel), emission of lines with E up ^ 30 000 
K is present only along the bow flanks, while lines with E„ p 
^ 30 000 K are observed in the whole shock, peaking at the 
bow head. Therefore, two main results emerge: (i) molec- 



1 The reduced images are available at 

http:/ / cdsweb. u-strasbg.fr 



ular hydrogen survives also where ionic emission is strong 
(see below) and, (ii) temperature gradients do exist along 
the shock, with the highest values reached at the bow head, 
where stronger excitation conditions are expected to occur. 
Atomic lines are listed in Table [5] and some examples of 
the observations are shown in Figure [1] middle and bottom 
panels. Plenty of [Fen] lines are detected (34 lines), emit- 
ted from levels with E up ^ 30 000 K. As for H2, two groups 
of lines are identified: those with E up S5 13 000 K, which 
come for the a 4 D level, are observed in the whole region, 
while those at higher excitation energy are emitted only at 
the bow head. In this same area, emission of hydrogen and 
helium recombination lines (8 and 3 lines, respectively, see 
Figure [TJ bottom panel) along with fine structure lines of 
[P 11], [Co 11], and [Tin] are detected. Other fine structure 
lines commonly observed in Herbig-Haro objects (e.g. [C 1] 
at 0.98/im, [Ni] at 1.04^m, [S11] at 1.03/xm, Nisini et al., 
2002) are not covered with SINFONI. In fact, the first wave- 
length in the J band is A=1.10/im, i.e. longer than that of 
other infrared spectrographs (e.g. ISAAC and SofI at ESO 
are able to observe wavelengths longer than 0.98/mi and 
0.86^m, respectively). 

To check for possible line variability, the intensities of 
atomic and molecular lines observed with SINFONI were 
compared with the ones observed by MC04 in July 2002 (i.e. 
four years before our SINFONI observations). Synthetic 
slits corresponding to slit 1 and 2 of MC04 (PA=32.4° and 
329.5°, width 0.6") were superimposed onto the SINFONI 
images and the flux was integrated over the same regions 
as in that paper. The main difference between SINFONI 
and ISAAC spectra is represented by the significant larger 
number of lines observed with SINFONI (see Figurcs[2]|4] 
and Figure 3 and Table 1 of MC04). In particular, it is no- 
ticeable the detection of many lines at high excitation, both 
atomic and molecular, that indicate excitation conditions 
stronger than those inferred in that paper (see Table[3]). 
A general increase in the intensities of lines observed with 
both the instruments is registered, but with differential be- 
haviour among molecular and ionic lines, which are respec- 
tively 4 and 1.5 times brighter than four years before. Such 
variability has been observed in proper motion studies, over 
time periods of a few years typical of radiative cooling times 
in HH objects (e.g. Caratti o Garatti et al., 2008). 

3. Diagnostics of physical parameters 

3.1. Fe analysis 

3.1.1. On [Fen] Einstein coefficients 

In Table[2] the [Fe 11] lines are listed grouping transitions 
originating from the same upper level. The intensity ratio of 
pairs of (optically thin) lines in each group is independent of 
the local physical conditions, being function only of atomic 
parameters (i.e. line frequencies and Einstein A coefficients 
for spontaneous emission). Hence, the observed intensity 
ratios among these lines can be efficiently used to measure 
the local extinction (e.g. Gredel, 1994). Unfortunalcty, the 
complexity of the energy level system of iron makes very 
difficult to accurately compute the A values, so that three 
distinct sets of these parameters, differing by more than 
30%, have been sofar listed in the literature: two were com- 
puted with different methods by Quinet, Le Dourneuf & 
Zeippen, (Q-SST, Q-HFR, 1996), and one is provided by 
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1E-18 3E-18 5E-18 1E-18 2E-18 3E-18 



Fig. 1. Selected spectral images from our HH99B data-cube. Intensities are given in colour scale. Offsets are from 
a 2 ooo = 19' l 02 m 05.4 ;s , <5 2 ooo=-36°54'39". A) H 2 :l-0 S(l) at 2.122^m. The locations of the knots labelled by D99 (Bl, B2 
and B3) and MC04 (BO) are indicated; B) H 2 : 2-lS(17) at 1.758pm; C) [Fen]: & 4 D 7/2 -& i F g/2 at 1.644pm. The black 
line delimits the area where [Fen] lines at S/N > 100 are detected and used to construct the plot of Figure^ D) [Fen]: 
a 4 P 3 / 2 -a 4 L> 7/2 at 1.749pm; E) H:Pa/3 at 1.282pm; F) [Pn]: 2 D 2 - 3 P 2 at 1.188/xm. White rectangles in panels A) and C) 
indicate the areas over which we have integrated the signal and extracted the spectra shown in FiguresHJIU 
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Fig. 2. J band spectrum extracted in knot B3 (black) and BO (red), over an area of 0.5 and 1.25 arcsec 2 , respectively. 



Nussbaumer & Storey (NS, 1988). A fourth list, based on 
the observation of P Cygni, was recently published by Smith 
& Hartigan (SH, 2006) who find A values 10%-40% higher 
than theoretical computations. The application of one set 
rather than another can lead to significant differences in 
the derivation of the local extinction. For example, the NS 
coefficients of the 1.644/xm and 1.257//m lines (both coming 



from the level & 4 D 7 / 2 with E„ p =11446 K and more com- 
monly observed) provide an extinction 2.7 mag higher than 
that derivable from the Q-SST coefficients. This marginally 
affects the NIR lines (the intensity grows by a factor of 2.8 
for a line at l^tm and a by factor of 1.4 for a line at 2/im), 
but becomes critical at optical wavelengths (the intensity 
grows by a factor of 33.6 at 0.5/^m). 
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Fig. 3. As in Figure [2] for the H band. 



The large number of [Fen] lines detected in HH99B, ob- 
served with a high S /N ratio in a remarkably large fraction 
of pixels, offers us the opportunity to compare theoretical 
predictions on the spontaneous emission rates with a sig- 
nificant sample of observational points. 

To this aim, we have plotted in Figure[5] the ratios 
I(1.257^m)/I(1.644 / um) vs. I(1.321/im)/I(1.644/^m), since 
these are observed at very high S/N (larger than 100, red 



filled squares). In the same Figure, green dashed curves rep- 
resent the Rieke & Lebofsky (1985) extinction law|3 applied 



2 The following arguments remain valid if other extinction 
laws are adopted (e.g. Cardelli, Clayton & Mathis, 1988), since 
these do not appreciably differ from the Rieke & Lebofsky law 
over the short wavelength range taken into consideration (from 
1.257/mi to 1.644 ^m). 
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Fig. 4. As in Figure [2] for the K band. 



to the intrinsic ratios expected for the four sets of A coeffi- 
cients (Q-SST, Q-HFR, NS, and SH). Squares along these 
'extinction curves' indicate Ay = 0, 5, 10 mag. 

First of all, we note that all the HH99B data lie defini- 
tively to the right of any of the plotted extinction curves. 
Since different Ay can move the points only along ex- 
tinction vectors, no theoretical intrinsic ratio is consistent 
with the observed points. This result was already pointed 



out in Nisini et al. (2005), who discussed how the extinc- 
tion along the knots of the HH1 jet determined from the 
1.321/im/1.644/im ratio is always smaller than that de- 
rived from the 1.257/im/1.644/im ratio, irrespective of the 
adopted theoretical set of A coefficients. HH99B data are 
also inconsistent with the P Cygni datum (green triangle in 
Figure [5]) , and consequently with the A coefficients extrap- 
olated from it. 
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Table 1. H2 lines observed in HH99B (to be continued). 



H2 lines 



Line id. 


^vac 

(H 


(K) 




Line id. 


(^m) 


(K) 




Line id. 


(urn) 


E U p 

(K) 














v= 


1 












l-OS(O) 1 


2.2233 


6471 


44 


1-0S(7) 


1.7480 


12818 


65 


1-0S(21) 3 


1.7195 


38136 


7 


1-0S(1) 


2.1218 


6951 


650 


1-0S(8) 


1.7147 


14221 


70 


1-0S(23) 


1.7801 


42122 


4 


1-0S(2) 


2.0338 


7584 


48 


1-0S(9) 


1.6877 


15723 


46 


1-0Q(1) 


2.4066 


6149 


44 


1-0S(3) 
1-0S(5) 2 


1.9576 


8365 


46 


1-0S(10) 


1.6665 


17312 


9 


1-0Q(2) 


2.4134 


6471 


30 


1.8358 


10342 


5 


1-0S(18) 3 


1.6586 


32136 


3 


1-0Q(3) 


2.4237 


6951 


50 


1-0S(6) 


1.7879 


11522 


30 


1-0S(19) 3 


1.6750 


34131 


8 


1-0Q(4) 


2.4375 


7584 


18 



v=2 



2-0S(0) 


1.1896 


12095 


6 


2-0Q(8) 
2-0Q(9) 5 


1.3020 


16881 


9 


2-lS(l) 


2 


2477 


12551 


44 


2-0S(l) 


1.1622 


12551 


48 


1.3188 


18108 


12 


2-lS(2) 


2 


1542 


13151 


22 


2-0S(2) 


1.1382 


13151 


7 


2-0Q(10) 


1.3375 


19436 


3 


2- IS (3) 


2 


0735 


13891 


40 


2-0S(3) 


1.1175 


13891 


6 


2-0Q(15) 


1.4648 


27267 


3 


2- IS (4) 


2 


0041 


14764 


21 


2-0S(4) 


1.0998 


14764 


6 


2-0Q(19) 3 


1.6156 


34447 


3 


2-lS(5) 


1 


9448 


15764 


30 


2-0Q(l) 


1.2383 


11790 


14 


2-00(2) 


1.2932 


11636 


7 


2-lS(8) 6 
2-lS(9) 7 


1 


8183 


19435 


7 


2-0Q(2) 4 


1.2419 


12095 


9 


2-00(3) 


1.3354 


11790 


16 


1 


7904 


20855 


22 


2-0Q(3) 


1.2473 


12551 


22 


2-00(6) 


1.4870 


13151 


5 


2-lS(10) 


1 


7688 


22356 


3 


2-0Q(4) 


1.2545 


13151 


7 


2-00(7) 


1.5464 


13891 


8 


2-lS(15) 8 


1 


7387 


30794 


6 


2-0Q(5) 


1.2636 


13891 


11 


2-00(9) 


1.6796 


15764 


5 


2-lS(17) 3 


1 


7587 


34446 


4 


2-0Q(6) 


1.2745 


14764 


7 


2-00(11) 


1.8349 


18108 


3 


2-lS(18) 3 


1 


7753 


36301 


3 


2-0Q(7) 


1.2873 


15764 


9 


2-lS(0) 


2.3556 


12095 


10 


2-lS(19) 3 ' 9 


1 


7962 


38166 


8 



v=3 



3-lS(0) 10 


1.2621 


17388 


9 


3-lS(9) 


1.1204 


25661 


22 


3-10(5) 


1.5220 


17819 


10 


3-lS(l) 


1.2330 


17819 


10 


3-lQ(l) 


1.3141 


17099 


8 


3-10(7) 14 


1.6453 


19087 


9 


3-lS(2) 


1.2077 


18387 


8 


3-lQ(2) 12 


1.3181 


17388 


12 


3-10(9) n 


1.7898 


25661 


10 


3-lS(3) 


1.1857 


19087 


9 


3-lQ(3) 


1.3240 


17819 


8 


3-2S(l) 


2.3864 


17819 


8 


3-lS(4) 


1.1672 


19913 


8 


3-lQ(5) 13 


1.3420 


19087 


9 


3-2S(2) 


2.2870 


18387 


9 


3-lS(5) 


1.1520 


20857 


10 


3-lQ(ll) 


1.4479 


25661 


8 


3-2S(3) 


2.2014 


19087 


25 


3-lS(6) 


1.1397 


21912 


5 


3-lQ(13) 


1.5024 


28558 


6 


3-2S(4) 


2.1280 


19913 


11 


3-lS(7) 


1.1304 


23071 


7 


3-lQ(15) 3 


1.5685 


31691 


6 


3-2S(5) 


2.0656 


20857 


17 


3-lS(8) 


1.1241 


24323 


5 


3-10(4) 


1.4677 


17388 


6 


3-2S(7) 


1.9692 


23071 


14 



To be confident of the reliability of our data, we have 
checked for possible blendings of the [Fe 11] lines with other 
lines or telluric features. In this respect, the 1.644/im line 
is close to both the H2 line 3-10(7) at 1.6453/im (see 
Table[T]and Figure[3]) and a telluric OH feature at 1.6442/im 
(Lidman & Cuby, 2000). This latter has been removed in 
the sky subtraction procedure, and its residuals have been 
estimated to affect the 1.644/im flux by less than 0.5%. The 
H2 line, coming from a low excitation level (E up ~ 19000 
K), is observed only in the receding parts of the bow, and 
thus it results spatially separated from the 1.644/Ltm emis- 
sion region. 

To check for other observational or data-reduction biases 
(e.g. unfavourable observational conditions, flat-fielding, in- 
tercalibration of lines lying in different bands), we searched 
in the literature for other observations of the consid- 
ered lines obtained with other instruments. To minimise 
the uncertainties, we have considered only line ratios ob- 
served with S/N > 30, that are shown with different 
colours/symbols in Figure Notably, all of them occupy 
the right side of the plot, in agreement with the HH99B 
points. This result, that reinforces the reliability of our 
observations, allows us to derive new A coefficients from 
our observations, provided that these are accurately cor- 



rected for the visual extinction value (measured indepen- 
dently from [Fen] lines). In this respect, two facts have 
to be noted: (i) although the sky area considered for this 
analysis is a few arcsec 2 (marked in black in Figure [1] mid- 
dle left panel), an extinction gradient of ~ 1 mag occurs 
in this zone, as evidenced by the scatter among the data 
points (red squares) of Figure [5J (ii) in the same area, we 
are able to obtain just a gross estimate of Ay (1.8±1.9 
mag) from the observed Pa/Br7 ratio (see Sect. 3.3). Both 
these circumstances prevent us from deriving an accurate 
measure of the A ratios. A rough estimate can be how- 
ever obtained de-reddening the average of the HH99B data 
(black cross) for Ay = 1.8 mag, which gives A1.321/A1.644 
= 0.38 and A1.257/A1.644 = 1.24. The main uncertainty on 
these values comes, more than from the error on the line 
fluxes, from the Ay estimate. Taking into account the Ay 
error of 1.9 mag, we can state that the theoretical points 
do belong to the segment of the extinction curve (along 
which they are constrained to move) starting at the point 
[0.33,1.02] (Ay = mag) and ending at the point [0.44,1.48] 
(Ay = 3.7 mag) . 

In conclusion, from a pure observational point of view, 
we can summarise as follows: (i) all the theoretically de- 
rived A values fail to reproduce the large majority of the 



8 

Table 1. H 2 lines observed in HH99B (continued). 



H2 lines 



Line id. 


^vac 


Eup S/N ma:c 


Line id. 




E U p S/N ma x 


Line id. 




E U p S f~N ma x 




w 


(K) 




w 


(K) 




(/irn) 


(K) 



u=4 



4-2S(0) 15 


1.3425 


22354 


9 


4-2S(8) 


1.1987 


28885 


7 


4-20(4) 18 


1.5635 


22354 


6 


4-2S(l) 


1.3116 


22760 


14 


4-2S(9) 


1.1958 


30141 


6 


4-20(5) 


1.6223 


22760 


7 


4-2S(2) 


1.2846 


23296 


21 


4-2Q(7) 


1.4592 


25625 


6 


4-20(6) 


1.6865 


23297 


4 


4-2S(3) le 


1.2615 


23956 


9 


4-2Q(9) 


1.4989 


27708 


7 


4-20 (7) 


1.7563 


23956 


7 


4-2S(4) 17 


1.2422 


24735 


9 


4-2Q(ll) 3 


1.5495 


30141 


3 


4-3S(3) 


2.3445 


23956 


6 


4-2S(5) 


1.2263 


25625 


12 


4-2Q(13) 3 


1.6123 


32857 


6 


4-3S(4) 


2.2667 


24735 


3 


4-2S(6) 


1.2139 


26618 


5 


4-2Q(15) 3 


1.6892 


35786 


6 










4-2S(7) 


1.2047 


27708 


10 


4-20(3) 


1.5099 


22081 


6 




















v= 


5 












5-3S(3) 


1.3472 


28500 


4 


5-3Q(l) 


1.4929 


26737 


3 


5-3Q(7) 3 ' 19 


1.5626 


30065 


6 


5-3S(4) 


1.3270 


29231 


3 


5-3Q(2) 


1.4980 


26994 


9 


5-3Q(ll) 3 
5-3Q(12) 3 


1.6673 


34291 


4 


5-3S(5) 


1.3107 3 


30066 


4 


5-3Q(4) 


1.5158 


27880 


5 


1.7021 


35529 


4 


5-3S(7) 


1.2894 3 


32017 


7 


5-3Q(6) 


1.5443 


29230 


3 


5-3Q(13) 3 


1.7412 


36821 


5 












v= 


6 












6-4S(l) 3 


1.5015 


31664 


3 


6-4Q(7) 3 


1.6829 


34175 


5 


6-40(4) 3 ' 21 


1.7965 


31306 


8 


6-4Q(3) 3 


1.6162 


31664 


5 


6-4Q(9) 3 - 20 


1.7369 


35992 


6 




















v= 


7 












7-5Q(5) 3 


1.7784 


36591 


3 



















Notes: "maximum signal-to-noise ratio in the line image. In case of blends the reported number refers to the sum of the blended 
lines (unless the emission comes from different zones of the bow). 



blends with [Fen] a 2 H 11/2 -a 2 G 



*9/2 

contaminated by atmospheric absorption 

detected in the whole bow 

blends with 4-2S(4) 

blends with 3-lQ(2) 

blends with HI 4-9 

blends with 3-10(9) 



8 blends with 6-4Q(9), HI4-10 

9 blends with 6-40(4) 

10 blends with 4-2S(3) 

11 blends with 2-lS(9) 

12 blends with 2-0Q(9) 

13 blends with 4-2S(0) 

14 blends with [Fell] a 4 D 7/2 



-a 4 F 



9/2 



15 blends with 3-lQ(5) 

16 blends with 3-lS(0) 

17 blends with 2-0Q(2) 

18 blends with 5-3Q(7) 

19 blends with 4-20(4) 

20 blends with 2-lS(15), HI 4-10 

21 blends with 2-lS(19) 



observed line ratios, irrespective of the extinction values; 
(ii) the best 'recipe' to derive a reliable extinction estimate 
from [Fen] lines is (at least when only the three consid- 
ered lines are detected) to use the NS coefficients for the 
1.321/j,m/1.644/im ratio and the Q-HFR coefficients for the 
1.257/im/1.644/imratio, that are 8% and 5% lower than our 
determinations; (Hi) dedicated observations of objects with 
well known visual extinction should be performed to derive 
the [Fe 11] Einstein A coefficients with sufficient accuracy. 

3.1.2. Extinction map 

Given the problems with the A coefficients outlined in the 
previous section, we have applied the following procedure 
to construct an extinction map across the HH99B bow from 
the observed [Fe 11] lines: to minimise the effect of the un- 
certainties we used a number of line ratios involving bright 
lines from four energy levels (i.e. a 4 D 7 / 2 , a 4 D 5 / 2 , a 4 D 3 / 2 , 
a 4 P 5 / 2 ) and distant in wavelength. With this set of ratios, 
and adopting, as a first attempt the NS coefficients, we 
determined the extinction in a very small region at the 
bow- head, where all the lines are detected with S/N ratio 



larger than 30. This value has then been used to calibrate 
the extinction map obtained from the 1.25/xm/1.644/im ra- 
tio, which is the only one detected well above the noise 
level (at least at 5<r) also in the bow flanks. Contours of 
the final map are shown in Figure[S] in a total area of ~ 
10 arcsec 2 variations of Ay up to 4 mag are recognised. The 
highest Ay values (4-5 mag) are found at the bow-head: 
if we correct the observed [Fe 11] lines for them, the emis- 
sion peak (see Figure [7]) moves of about 0.6 arcsec towards 
north-east. Along the flanks Ay is generally lower (up to 
2-3 mag) . Thus, the progressive fading in these zones of the 
Fen emission can not be ascribed to an increasing extinc- 
tion, but rather reflects low excitation conditions and/or 
low abundance of the gas-phase iron. 

As described above, the main uncertainty on the extinction 
map arises from the adopted set of the A values. We have 
thus re-derived the same map from the 1.257/im/1.644/im 
ratio, now adopting our Einstein coefficients ratio of 1.24. 
We find that the largest difference between the two maps 
occurs at the emission peak, where it is of ~ 0.6 mag. This 
implies a marginal increase in intrinsic line intensities (for 
example I(1.644^m) increases by 10%) and does not crit- 
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Table 2. Ionic lines observed in HH99B. 



Ionic lines 



Line id. 




E^p S/N ma3; 


Line id. 


^vac 


E U p S f~Nmax 




(H 


(K) 




w 


(K) 



[Fell] lines'' 



a 4 D 7/2 


-a°D 9/2 


1 


2570 


11446 


700 


a 4 D 1/2 -a°D3 /2 


1 


2525 


12729 


20 


a 4 D 7/2 - 


a 6 D 7 /2 


1 


3209 


11446 


130 


a 4 Di /2 -a 6 Di /2 


1 


2707 


12729 


24 


a 4 D 7/2 - 


a 4 F 5/2 


1 


9541 


11446 


7 


a 4 D 1/2 -a 4 F 5/2 


1 


6642 


12729 


18 


a 4 D 7/2 - 


a 4 F 9/2 1 


1 


6440 


11446 


720 


a 4 D 1/2 -a 4 F 3/2 


1 


7454 


12729 


21 


a 4 D 7/2 - 


a 4 F 7/2 


1 


8099 


11446 


99 


a 4 P 5 /2-a 4 D 5 / 2 


1 


9675 


19387 


20 


a 4 D 5/2 


-a^g/a 


1 


1916 


12074 


7 


a 4 P 5/2 -a 4 D 1/2 2 


2 


1609 


19387 


5 


a 4 D 5/2 - 


a 6 D 7/2 


1 


2489 


12074 


17 


a 4 P 3 / 2 -a 4 D 7/2 


1 


7489 


19673 


9 


a 4 D 5/2 - 


a 6 D 5 /2 


1 


2946 


12074 


70 


a 4 Pi/ 2 -a 4 D 5/2 


1 


8139 


20006 


62 


a 4 D 5/2 - 


a^s/a 


1 


3281 


12074 


42 


a G 9 / 2 — a D 7 / 2 


1 


2675 


22797 


7 


a 4 D 5/2 - 


a 4 F 9/2 


1 


5339 


12074 


30 


a ijr 7 / 2 — a D7/2 


1 


1885 


23552 


16 


a 4 D 5/2 - 


a 4 F 7/2 


1 


6773 


12074 


157 


a 2 P 3 /2-a 4 P 5/2 


2 


0466 


26417 


4 


a 4 D 5/2 - 


a 4 F 5/2 


1 


8005 


12074 


35 


a 2 P3/2-a 4 P 3 / 2 


2 


1334 


26417 


5 


a 4 D 3 / 2 


-efD 3 /2 


1 


2791 


12489 


41 


a 2 P 3 /2-a 4 Pi/ 2 


2 


2442 


26417 


4 


a 4 D 3/2 - 


at)i/a 


1 


2981 


12489 


23 


a 2 H 11 / 2 — a 2 G 9 / 2 4 


2 


2244 


29265 


21 


a 4 D 3/2 - 


a 4 F 7/2 


1 


5999 


12489 


46 


a 2 G 9 / 2 


2 


0157 


29934 


10 


a 4 D 3/2 - 


a 4 F 5/2 


1 


7116 


12489 


13 


a 2 Hg/2— a 2 G 7 / 2 


2 


2541 


29934 


5 


a 4 D 3/2 - 


a 4 F 3/2 


1 


7976 


12489 


36 


b 4 Pi /2 -a 4 P 3/2 


1 


1446 


32242 


4 



H lines 



3-5 (Pa/3) 


1.2822 


151492 


21 


4-11 


1.6811 


156499 


6 


4-14 


1.5884 


156999 


3 


4-10 5 


1.7367 


156226 


6 


4-13 


1.6114 


156870 


4 


4-9 6 


1.8179 


155855 


6 


4-12 


1.6412 


156708 


6 


4-7 (Bry) 


2.1661 


154583 


12 



Other lines 



He^Pi - l S 


2.0587 


246226 


3 


[Tin 


a 2 F 5/2 -a 4 F 3/2 10 


2.1605 


6652 


5 


Hei 3 L>- 3 P 7 


2.0607 


282101 


8 


[Tin 


a 2 F 7/2 -a 4 F 9/2 


2.2201 


7040 


5 


Hei 3 P - 3 D 7 


1.9522 


289992 


5 


[Till 


a 2 £> 5 /2-a 4 F 5 /2 


1.1560 


12570 


3 


[P n] 2 D 2 - 3 P 1 


1.1471 


12764 


5 


[Till 


a 4 P 3/2 -b 4 F 7/2 


1.2036 


13506 


3 


[Pn] 2 D 2 - 3 P 2 8 


1.1886 


12764 


16 


[Till 


a 2 P 3 /2-b 4 P 3/2 


1.1028 


14340 


3 


[Co n] b 3 F 4 -a, 5 F 5 9 


1.5474 


14119 


5 


[Till 


a 2 #n/ 2 -a 2 F /2 


1.2695 


17729 


5 



Notes: "maximum signal-to-noise ratio in the line image. In case of blends the reported number refers to the sum of the blended 
lines (unless the emission comes from different zones of the bow). b Lines coming from the same upper level are grouped, and the 
first term of each group is evidenced with bold-face characters. 

1 blends with 3-10(7) 

2 blends with [Ti n]a 2 P 5/2 -a 4 F 3/2 (fundamental transition) 

3 blends with [P ll] 1 D 2 - 3 P 2 

4 blends with 1-0S(0) 

5 blends with 2-lS(15), 6-4Q(9) 



6 blends with 2-lS(8) 

7 multiplet 

8 blends with [Fell] a 2 G 7/2 -a 4 D 7/2 

9 tentative identification 

10 blends with [Fell] a 4 P 5/2 -a 4 D 1/2 



ically affects the derivation of the physical parameters of 
the atomic gas (see next section). 

3.1.3. Electron density map and temperature 

To derive the electron density along the bow struc- 
ture, we selected seven intensity ratios (i.e. Ji, 533/^1.6447 

-^1.60o/^l. 644, II. 677 /II. 644; A. 664 /A. 644, ^1.27l/A.257, 

h. 279/^1.257, J1.328/A.257)) involving lines close in wave- 
length (their dfferential extinction is negligeable) and 
coming from levels with different critical densities (from 
~ 8 10 2 to 3 10 5 cm -3 ) and similar excitation energy 
(E up - 11000-12 000 K), so that the dependence on the 



temperature is very weak. All these line ratios have been 
simultaneously fitted with a NLTE code which solves the 
equations of the statistical equilibrium for the first 16 fine 
structure levels of [Fe n] . Spontaneuos rates are taken from 
NS, while energy levels and rates for electron collisions 
are adopted from Pradhan & Zhang (1993). Assuming 
T e = 10 000 K, we have constructed the electron density 
map shown in Figure[8l n e is of the order of 2-4 10 3 cm~ 3 , 
with a peak up to 6 10 3 cm -3 at the bow heacEl These 
values are in the range commonly found in HH objects 

3 As we will see in the following, we estimate in this part of 
the bow T e ~ 16000 K. At this temperature the ratios involved 
in the electron density estimate can differ by less than 15% with 
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- ■ HH99 - SINFONI 

" ° HH111 - SofI ©H 
" * HH240 - Sofl/ISAAC 




0.1 0.2 0.3 0.4 



I(1.321/xm)/I(1.644/xm) 



Fig. 5. I(1.257/im)/I(1.644/im) vs. I(1.321/.tm)/I(1.644/im) ratios measured in different objects (depicted with different 
symbols/colors). The HH99B-SINFONI data (red squares) have been computed in pixels where the S/N of each of the 
three lines is larger than 100, while other data are literature observations where S/N > 30. Intrinsic line ratios predicted 
theoretically (Q-SST = Quinet et al. 1996 - Superstructure; Q-HFR = Quinet et al., 1996 - Relativistic Hartree-Fock; 
NS = Nussbaumer & Storey, 1988), along with the observational point (SH) by Smith & Hartigan (2006), are labelled. 
Green dashed curves represent the extinction law by Rieke & Lebofsky (1985), starting from different theoretical points; 
open squares refer to Ay =0,5,10 mag. The same extinction law (in red) has been applied to the Ay = mag point 
derived from SINFONI data. This latter has been derived by applying to the average of the HH99B points (black cross) 
a visual extinction of 1.8 mag, as estimated from the Pa/Br7 ratio (see text). References: HH99B - SINFONI: this work; 
HH111-, HH240-, HH120-Sofl: Nisini et al., 2002; HH240- ISAAC: Calzoletti et al., 2008; HH99B - ISAAC : M c Coey et 
al., 2004; HH1 - SofI: Nisini et al., 2005; HH54 - ISAAC: Giannini et al., 2007; HH34 - SofI: Podio et al., 2006; HH135 - 
SofI: Gredel, 2006; Orion bar - SofI: Walmsley et al., 2000; P Cygni - Spex: Smith & Hartigan, 2006. 



from embedded jets (e.g. Nisini et al., 2005, Podio et al., 
2006). 

In a restricted area at the bow head of about 1 arcsec 2 
(see Figured! middle right panel), we have detected 13 lines 
at high excitation (E up between 20 000 and 30 000 K), 
which are suitable for evaluating the local electronic tem- 
perature. Of these, just four lines coming from the term a P 
(i.e. 1.749/xm, 1.814/mi, 1.967/^m, 2.161^im) can be mod- 
elled, since for the remaining nine lines the collisional rates 
are unknown. 

To that aim, we have enlarged our 16 level code by in- 
cluding further three fine structure levels for which the colli- 
sional coefficients are reported by Zhang & Pradhan (1995). 
Notably, the excitation energy of level #19 is around 32 000 
K, well above than that of level #16 (less than 20 000 K): 
hence the temperature range that can be probed with the 



19 level code is sensitively enlarged. Having fixed extinction 
and electron density from the maps of Figures [5] and [5J we 
fitted the de-reddened ratios with the 1.257/im line, inte- 
grated over the area where the a 4 P lines are detected at S /N 
>5. Results are plotted in Figured where the observed ra- 
tios are compared with the predictions of both the 16 and 
19 level codes and for temperatures from 10 000 to 20 000 
K. First, we note that while at T e ^10 000 K the inclusion 
of three further levels does not change the results of the 
16 level code, strong differences emerge at higher tempera- 
tures (e.g. the ratio 1.257/1.749 decreases by about 70% at 
T e =20 000 K). Second, ratios with the 1.814/zm, 1.967/zm, 
and 2.161/xm lines well agree with T e ~ 16 000 - 17 000 K 
(the latter line has been de-blended from the fundamental 
line of [Tin], using a NLTE model for this species, Garcia- 
Lopez et al., 2008). 



respect to those computed at T= 10000 K. This implies a very 
marginal increase of n e . 



Finally, we note that the ratio with the 1.749/im line im- 
plies T e ~8 000 K, that we consider not reliable because in 
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-4 -2 2 4 

offset (arcsec) 

Fig. 6. Extinction map obtained from [Fen] lines where contours from Ay=l mag to Ay =4 mag are shown. The NS 
Einstein A coefficients have been adopted. 




Fig. 7. . De- reddened contours of 1.257/mi intensity, overlaied with the image acquired in the same line. 



the same spatial region examined here also bright hydrogen 
and helium recombination lines are observed (see Figure[T]) . 
For this line, however, neither evident discrepancies in the 
different computations of the Einstein coefficients (all the 
available lists give similar values, see Sect. 3. 1.1), nor obser- 
vational biases (e.g. extinction, blending with other lines) 
are able to explain the disagreement with the other ratios. 

3.1.4. [Fen] abundance 

The gas-phase Fe abundance x(Fe ) is a measure of the 
shock efficiency in disrupting the cores of the dust grains 
where iron is locked in quiescent conditions (e.g. May et al., 



2000). Estimates of x(Fe ) in shock environments so far have 
given sparse results, from values close to solar abundance 
(e.g. Beck-Winchatz, Bohm & Noriega-Crespo, 1996), up to 
intermediate (Nisini et al., 2002) and very high depletion 
factors (Mouri & Taniguchi, 2000, Nisini et al., 2005). A 
powerful way to estimate the percentage of gas-phase iron 
(Sfc) based on [Fen]/[Pn] line ratios has been proposed by 
Oliva et al. (2001). Since phosphorus and iron have simi- 
lar ionisation potentials and radiative recombination coef- 
ficients, they are expected to be in the first ionised state 
in comparable percentages; moreover, the near-IR lines of 
Fen and Pll have similar excitation energies and critical 
densities and therefore are excited in similar physical con- 





Fig. 9. [h.257 / I\]obs/[h.257 / I\]mod plotted vs. the electronic temperature for four [Fen] lines. The results for the Fe 
NLTE model at 16 and 19 levels are shown for comparison. 



ditions. Hence, [Fe n]/[P n] line ratios are good indicators of 
the relative abundance of the two species, and more specifi- 
cally, because phosphorus is a non-refractory species, of the 
degree of iron depletion. In HH99B two [P n] lines are de- 



tected, at 1.1471 and 1.1885 /im, this latter barely blended 
with a [Fe n] line (see Table [lj . Oliva et al. derive: 



x(Fe) /([Fen] 1.257) _ /([Fell] 1.257) 

x(P) X 7([Pn]1.188) ~ 7([Pii]1.147) 
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offset (arcsec) 

Fig. 10. Map of the percentage of iron in gas-phase with overlaied the intensity contours of the [Fen] 1.257/im line. 



This equation, as stated by the authors, is accurate to 
within a factor of 2 for all temperatures and densities ex- 
pected within the shock. By assuming a solar Fe /P abun- 
dance ratio of ~ 120 (Asplund, Grevesse & Sauval, 2005) 
we have derived a map of the percentage of gas-phase iron 
(see Figure fT0|) . A strong decrease in the percentage of gas- 
phase iron occurs from the bow-head (70%) towards the 
zones behind (up to 20%), with an average uncertainty of 
about 15%, estimated from the propagation of the errors in 
both the two considered images. Notably, theoretical pre- 
dictions for the degree of iron depletion as a function of the 
shock velocity (Jones, 2000) imply that this latter should 
exceed 100 km s _1 for 5fc ^ 0.20, a condition verified in 
our case even in the bow flanks. 



3.2. H2 analysis 
3.2.1. Extinction map 

Looking at the line images of FigureQ]it is evident that the 
bulk of the H2 emission comes from the bow flanks. Since 
these parts of the bow are not covered in the extinction map 
obtained with [Fell] lines (see FigureH]), this latter cannot 
be used to de-redden the H2 emission. Hence, as a first step 
of the H2 line analysis, we have derived a new extinction 
map. Out of a number of ratios of lines coming from the 
same upper level, we considered only the three ratios (1- 
0S(1)/1-0Q(3), 2-0S(l)/2-lS(l), 2-0Q(3)/2-lS(l)) of lines 
detected at S/N per pixel larger than 5 over the whole emis- 
sion region. Since the one at the largest S/N (1-0S(1)/1- 
0Q(3)), suffers from the poor atmospheric transmission at 
2.42/^m of the 1-0Q(3), we have used the other two ratios 



as calibrators and the 1-0S(1)/1-0Q(3) ratio to probe the 
differential extinction along the shocked region. The Rieke 
& Lebofsky (1985) extinction law was adopted. The final 
map is shown in Figure fTTl the covered zone somewhat com- 
plements the Ay map constructed from [Fe 11] lines, with a 
partial spatial overlap in the areas corresponding to knots 
B3 and Bl of D99. The Ay values typically range from 1 to 
4 mag, in substantial agreement with those inferred from 
[Fe 11] emission: therefore this result demonstrates that the 
reddening within the shock does not change significantly 
and that the Ay variations mainly arise from the fore- 
ground material inside the cloud. Our result contrasts with 
the general trend observed both in other Herbig-Haro ob- 
jects (e.g. Nisini et al., 2002, Giannini et al., 2004) and 
in HH99B itself (MC04), where the extinction computed 
from [Fe 11] lines is systematically higher than that probed 
with H2 lines. This may stem because in these cases the 
extinction was computed from the 1.257/zm/1.644/zm ratio, 
for which the A coefficients of NS have been adopted (see 
Sect.3.1.1). 



3.2.2. Temperature map 

The temperature of the molecular gas can be obtained 
from a Boltzmann diagram (e.g. Gredel, 1994), plotting 
ln(N V: j I 'gj) against E Vi j/k. Here N Vt j (cm -2 ) is the col- 
umn density of the level (v, J),E v j (K) its excitation en- 
ergy and g,j = (2 J + 1)(2J + 1) the statistical weigth (here 
we assume that 1=0,1 for para- and ortho-H2). If the gas is 
thermalised at a single temperature, the data points align 
onto a straight line, whose slope gives the gas tempera- 
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Fig. 11. Extinction map as derived from H2 line ratios. Contours from Ay = l to 4 mag are shown. 



ture. In this diagram, points that refer to transitions from 
the same upper level do overlay each other, once corrected 
for extinction effects. We applied this method to all the 
pixels of the H2 images, and the resulting map is shown in 
Figure [T2J where contours of temperature are given in units 
of 10 3 K. Three results can be taken from this map: (i) a 
temperature gradient from ~ 2 000 K up to 6000 K occurs 
from the receding parts of the shock towards the head. We 
underline that this gradient can be traced because of the 
very large number of H2 lines detected, which cover the 
Boltzmann diagram up to excitation energies of ~ 38 000 
K, therefore sensitively enlarging the dynamical range of 
temperatures typically probed with H2 near-infrared lines; 
(ii) two different behaviours in the Boltzmann diagram oc- 
cur between the bow-head and the flanks (see the insets 
in Figure [12]): H 2 appears fully thermalised (at T~ 5000 
K) at the bow-head (hence here the contours give directly 
the gas temperature), while a curvature exists among the 
points in the southern flank diagram, so that at least two 
temperature components can be traced. In these parts of 
the bow the contours indicate the average of these temper- 
atures. The outlined behaviour can be generalised to the 
whole bow structure and likely reflects, from one side, that 
fluorescence can be discarded as a possible excitation mech- 
anism, since it implies strong departures from thermalisa- 
tion (Black & van Dishoeck, 1997), and, from the other side, 
that different shock mechanisms are at work. In fact, the 
response of the level populations to the shock parameters 
can be seen in the Boltzmann diagram. Roughly speaking, 
an enhancement of the shock velocity increases the rate of 



collisions that vibrationally excite the H2 molecules and 
favours thermodynamical equilibrium: thus the thermalisa- 
tion observed along the head would testifies the presence of 
a fast shock. If this were the C-type shock should be 

favoured, otherwise H2 could not survive (at high temper- 
atures) against dissociation (e.g. Le Bourlot et al., 2002, 
Flower et al., 2003). The alternative possibility that the 
observed H2 emission arises from re-forming H2 onto dust 
grains can be reasonably discarded. In fact, the timescale 
of this process is of the order of 10 17 /n s -1 (Hollcnbach & 
McKee, 1979), being n the density of the gas after the com- 
pression due to the shock passage. For reasonable values of 
this parameter (10 4 -10 7 cm -3 ), and assuming a velocity of 
the downstream gas of the order of 10 km s" 1 , the distance 
covered by the gas at the bow-head should be between 10 16 
and 10 19 cm, definitively longer than the linear dimension 
of knot B0, i.e. ~ 5 10 15 cm; (Hi) from the same Boltmann 
diagrams a marked difference in the column density of the 
warm H2 (from the intercept of the straight line) emerges. 
For example, from the insets of Figure[T2l we derive N(H2) 
= 2 10 17 cm -2 and 3 10 16 cm -2 at the flank and at the 
bow-head, respectively. If the same shock lengths are as- 
sumed for these regions, it derives that about 85% of the 
molecular gas must get destroyed in the most excited part 
of the shock (although H 2 emission is still observed). 

3.3. H analysis 

At the bow-head we observe hydrogen recombination lines 
of the Brackett series along with the Pa/3 line. The S/N 
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E yJ (10= K) 

Fig. 12. Temperature map (in 10 3 K) as derived from H2 line ratios. As an example, we show the rotational diagrams 
in two points of the bow: while at the bow-head H2 appears thermalised at (T~ 5000 K), in the southern flank two 
temperature components co-exist, with an average temperature (corresponding to the point shown in the map) of ~ 2800 
K. 



ratio of these lines, but Br7 and Pa/3, is so low, however, 
that detailed modelling is prohibitive. Thus, we have used 
just the Pa/3/Br7 ratio to derive an independent estimate 
of the extinction. Assuming case B recombination (Storey 
& Hummer, 1995), we obtain Ay = 1.8±1.9 mag. 
More interesting parameters, i.e. the hydrogen fractional 
ionisation, x e , and the hydrogen post-shock density, n# 
= n e /x e , are obtainable from the observed intensity ra- 
tio [Feii]1.257/Pa/3. As described in detail in Nisini et 
al.(2002), under the assumption that iron is singly ionised, 
such a ratio can be expressed as: 

x rv m\ /[Fell]1.257\ _1 e[ Fe n]i.257 
x e = S Fe Fe/H — — 2 

u \ Pa/3 J epa/3 

with 5 Fe the gas-phase iron fraction with respect to the 
solar Fe abundance, (Fe /H) Q and £[foii]i.257 an d epa/3 (hi 
erg cm 3 s _1 ) the emissivities of the two lines, taking ep a f3 
from Storey & Hummer (1995). The above quantity was 
computed for the physical conditions derived at the bow- 
head, i.e. T=16 000 K, n e =6 10 3 cm" 3 , S Fe ~ 0.7 and 
having taken (Fe/H) Q = 2.8 10~ 5 (Asplund, Grevesse & 
Sauval 2005). This leads to x e ~ 0.4-0.5 and n# between 0.8 
and 1.4 10 4 cm -3 . In the receding parts of the shocks, the 
fractional ionisation cannot be directly computed, since we 
have no estimates for the electron temperature. However, 
under the reasonable assumption that T e ^ 10 000 K a 
sharp decrease of x e is expected: for example, for the re- 
gions where 5i?e=0.3-0.4, we find x e ~ 0.2-0.3. These es- 
timates are in agreement with those inferred along other 



Herbig-Haro objects through optical line diagnostics (e.g. 
Bacciotti & Eisloffel, 1996, Hartigan & Morse, 2007) and 
with those inferred in a number of jets by combining optical 
and infrared observations (Nisini et al., 2005, Podio et al., 
2006). 

4. Kinematical properties 

4.1. H 2 

In this section we intend to characterise the kinematical pa- 
rameters of the shock(s) in HH99B. This topic was already 
discussed by D99, and we compare here our results with 
those found in that work. The higher spectral resolution 
(R~ 15 kins -1 ) of D99, obtained with echelle spectroscopy, 
has revealed that the peak velocity of the H2 2.122/im line 
moves progressively from slightly blue-shifted values near 
the shock front towards red-shifted values in the flanks, and 
has been interpreted in the framework of a receding bow- 
shock oriented with respect to the line of sight of about 45° . 

Nominally, the spectral resolution of our K band obser- 
vations would not permit us to reveal variations of the order 
of those measured by D99. This limitation, however, is par- 
tially compensated by the very high S/N ratio at which 
we detect the 2.122/j,m line. A trend on both the peak ve- 
locity and on the profile shape can be followed along the 
bow structure, although we cannot give precise numerical 
estimates on the line parameters (v pea k, FWHM). Our re- 
sults are presented in FigurefTSI where the contours of the 
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Fig. 13. Local standard of rest (LSR) velocity map of the 2.122/im line peak, with superimposed line intensity contours. 
We indicate with a white line the diameter 2Kdis of the last cap beyond which the bulk of H2 is dissociated (sec Sect. 4. 2). 
Insets show the line profile observed at the bow head (top left), in the southern flank (bottom left) and at the bow centre 
(bottom right). Blue and red asymmetries are visible at the bow head and along the flanks, while the line is symmetric 
toward the centre. The instrumental profile, measured on OH atmospheric lines, is shown for comparison (red dashed 
line) . 



2.122/xm line intensity (de-reddened) have been superposed 
on the Vp e ak map. 

Overall, our results confirm those of D99: the line pro- 
file presents a blue-shifted component towards the shock 
front at the bow head (BO). The opposite occurs along 
the two flanks and especially along the edge of the Bl 
flank (not covered by the echelle spectra in D99), where 
the line peak is shifted of ~ +15 km s _1 with respect to 
the line profile at the head. Analogous to the spectra of 
D99, the 2.122/im profile does not show double peaked com- 
ponents, as generally expected for a parabolic bow struc- 
ture (Schultz, Burton & Brand, 2005) though it does be- 
come wider near the centre of the bow, where the opposite 
sides are seen in projection. Here the observed FHWM„(, S 
is 85-105 km s _1 , that, deconvolved with the instrumental 
profile, measured on atmospheric OH lines, roughly gives 
an intrinsic line width of ~ 20-40 km s _1 . The agreement 
with previous observations is also mantained along the bow 
flanks, where the profile width becomes narrower decreas- 
ing toward the spectral resolution limit on the intrisic width 
of ~ 20 km s _1 . As for the line peak, a sudden increase of 
FHWM b s is registered at the edge of the southern flank, 
where we measure up to ~ 115 km s _1 , i.e. an intrinsic 
width of ~ 70 km s . This is close to the maximum shock 
velocity (~ 80 km s _1 ) at which H2 can survive against dis- 
sociation, predicted by the C-shock model by Le Bourlot et 



al. (2002). This last topic will be further discussed in the 
next section. 

4.2. [Feu] 

We have performed the kinematical analysis of the ionic gas 
component on the two brightest [Fen] lines at 1.257 and 
1.644/im. Both appear resolved in velocity and give similar 
results for the line profile shape (which is Gaussian across 
the whole bow), the intrinsic line width, the peak position 
and the FWZI. This last result gives a direct measure of 
the shock speed at the bow apex (Hartigan, Raymond & 
Hartmann, 1987): we obtain Vbow = H5±5 km s , which 
agrees well with the prediction by D99 (80-120 km s _1 ) 
derived from the overall shape of the bow and the charac- 
teristics of the H2 2.122/im line profile. 
In Figure[T4lwe show the radial velocity map of the 1.644/im 
line: this appears red-shifted over the whole bow structure, 
the maximum shift occurring towards the image centre. We 
intepret this behaviour as a geometrical effect due to the in- 
clination of the bow with respect to the line of sight. Indeed, 
if the bow is observed at a certain angle a ^ 0°, 180°, the 
peak of the radial velocity component is seen offset from the 
bow apex. Following the procedure described in Appendix 
A, we are able to express the projected distance (D') be- 
tween the emission and the radial velocity peaks as a func- 
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Fig. 14. Left panel: local standard of rest (LSR) velocity map of the 1.644/im line peak, with superimposed line intensity 
contours. The line is red-shifted along the bow-structure, with a peak at about +70 km s _1 . The projection D' along 
the line of sight of the distance between the line emission peak and the maximum radial velocity is shown, together with 
the diameter 2R ( iis (see caption of Figure [T3"]) . Right panel: schematic view of the parabolic geometry assumed for the 
bow-shock. The apex coincides with the origin of the coordinate system. R=^/.x 2 + y 2 , is the radius at the distance z 
along the bow axis, <f> is the azimuthal angle around the z-axis, 9 is the angle between the bow direction and the tangent 
to the parabolic surface (at each the shock velocity at the apex, v&, is de-composed in V|| and vj_) and a is the angle, 
lying on the x-z plane, between the z-axis and the line of sight. D is the distance between the line emission peak and the 
maximum radial velocity and D' its projection onto the sky plane. 



tion of a and the H2 breakdown velocity (vdis), once the 
relationship between this latter and the parameter b regu- 
lating the bow (parabolic) shape is defined. The result is 
depicted in Figure [15] interestingly, we found a values close 
to that inferred by D99 for a range of Vdis between 70 and 
90 km s _1 . These values can be considered as an indirect 
measurement of a parameter whose theoretical predictions 
have been widely discussed over the last decades: a first 
value of 24 km s _1 was inferred by Kwan (1977), subse- 
quently Draine, Roberge & Dalgarno (1983) found Vdi S =50 
km s _1 , a result confirmed by Smith & Brand (1990) and 
Smith (1991). More recently, Le Bourlot et al. (2002) have 
shown that Vdis can increase up to 80 km s _1 for values 
of the pre-shock density around 10 3 cm -3 . Our measure- 
ment is therefore in agreement with this latter prediction. 
Moreover, v^is ~ 80 km s _1 is consistent with the FWHM 
of the 2.122/im line profile measured in the southern flank 
(Sect. 4.1), and reinforces the hypothesis that a fast, contin- 
uous shock is responsible for the excitation of the molecular 
gas at the apex of the bow (Sect. 3. 2. 2). 



5. Concluding remarks 

We have presented bi-dimensional, deep near-infrared spec- 
tral images of the bow shock HH99B. These have allowed 
us, for the first time, to accurately derive the physical pa- 
rameters of both the molecular and ionic gas components 
(summarised in Tabled where our results are compared 
with those derived in previous works), and, at the same 
time, to characterise the geometry and the kinematical 



properties of the flow. The main results of this study are 
the following: 



- More than 170 emission lines have been detected, mainly 
ro- vibrational H2 and [Fe 11] lines, many of them never 
observed before in an Herbig-Haro object. In addition, 
transitions of hydrogen and helium recombination and 
fine structure lines of [Pn], [Tin], and possibly [Coil] 
have been observed. 

- A clear bow-shape morphology emerges from the line in- 
tensity maps. As shown in Figures [2EJ [Fe 11] and other 
ionic emission peaks definitely at the bow-head (B0), 
being strong also in the knot B2 immediately behind. In 
contrast, H2 emission delineates the bow flanks, peaking 
in the knots Bl and B3. Notably, the H2 lines with the 
highest excitation energy (E up il 30 000 K) show a dif- 
ferent morphology, being strong towards the bow-head. 
This implies that H2 still survives in this zone, in spite 
of the significant temperature enhancement there. 

- Extinction maps have been derived from the analysis of 
both [Fe 11] and H2 lines. These give similar results, with 
Ay ranging between 1 and 4 mag. 

- A detailed electron density map has been obtained in 
the framework of NLTE approximation for [Fen] line 
emission. This remains almost constant in the [Fen] 
emission zone, peaking towards the bow-head. From the 
same emission, we are able to probe a variation of the 
electron temperature, which falls from ~ 16 000 K at 
the apex to less than 10 000 K in the receding parts of 
the bow. 
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Fig. 15. Projection over the sky plane of the distance between the [Feii]1.644/im emission peak and the maximum 
radial velocity, plotted against the bow inclination angle a. The hatched area defines the uncertainty on D' , estimated 
as the maximum and minimum distance between the 1.644//m line peak and the points where the maximum shift in 
velocity is registered (those coloured in red in Figure[T4l left panel). The curves refer to different assumptions for the H2 
breakdown speed, Vdisi having measured the shock velocity at the bow apex, Vbow an d the radius Hdis (see Appendix A). 
The observations are matched for Vdis between 70 and 90 km s -1 . 

Table 3. Physical parameters estimated in HH99B. 



Parameter This work MC04 D99 



Ay (Fen) (mag) 


1-4 


4-10 




Ay(H 2 ) (mag) 


2-4 


<4 


<4.6 


Ay(H) (mag) 


~2 






n e (cm 3 ) 


3-6 10 3 






Te (K) 


< 17 000 






Th 2 (K) 


2 500 - 6 000 


2 000 - 4 000 




Nh 2 (cm 2 ) 


3 10 16 - 2 10 17 






5p e il 


0.2-0.7 


0.25 




x e 


0.2-0.5 






riH (cm -3 ) 


~10 4 (post-shock) 


10 4 (pre-shock) 




v bow (km s -1 ) 


110-120 ([Fen]) 


50 ([Fen]) 


80-120 (Ha) 


inclination angle (°) 


40-60 




45 



- An iron depletion degree not higher than 30% has been 
inferred at the bow apex, which testifies in favour of a 
J-type shock as the main excitation mechanism in this 
part of the bow. In this same zone, we infer a fractional 
ionisation of ~ 0.6 and a post-shock density of ~ 10 4 
cm -3 . 

- Analysis of H2 line emission has allowed us to probe the 
molecular temperature variation. We find that at the 
bow apex thermalisation has been reached at T^ 6 000 
K, likely due to the presence of a fast, non-dissociative 
shock. On the contrary, along the flanks different tem- 
perature components are simultaneously present. A de- 
crease of about a factor of ten is registered in the H2 
column density going from knots B1/B3 toward knot 
B0. Both these circumstances are not accounted for by 
models which have attempted to interpret the H2 emis- 
sion on HH99B on the basis of fewer lines. Therefore, 



the conclusions of such models should be tested in the 
light of this new piece of information. 

- From the brigthest [Fen] and H2 lines we have been 
able to probe the kinematical properties (e.g. shock ve- 
locity) of the shocked gas. In particular, we confirm the 
result by D99 according to which HH99 is a red-shifted, 
receding bow. 

- The radial velocity map of [Fe 11] emission has been in- 
terpreted in the framework of the bow geometry. From 
this map we have consistently inferred the bow inclina- 
tion angle and defined a range of 70-90 km s -1 for the 
H2 breakdown velocity. These values are in agreement 
with the prediction of Le Bourlot et al.(2002) of Vdis up 
to 80 km s -1 . We propose our method as a valuable tool 
to derive the jet inclination angle (if larger than 10-20°) 
in cases where proper motion is unknown. 

- The kinematical parameters of the [Fe 11] emission esti- 
mated in this work do not confirm the model predictions 
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by MC04. In particular, the argument that the [Fen] 
lines originate in a pure J-type shock with v s hock ~ 50 
km s _1 contrasts with our measure of v s hock ~ 110-120 
km s _1 . Thus, the modelling of the [Fen] emission may 
have to be revised. 
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Appendix A: Derivation of inclination angle and Ho 
breakdown velocity 

With reference to Figure[T4l right panel, we assume that 
the bow has a parabolic geometry (z — R 2 /2b), where the 
apex coincides with the origin of the coordinate system. 
R= y/x 2 + y 2 , is the radius at the distance z along the bow 
axis, <t> is the azimuthal angle around the z-axis, 9 is the 



angle between the bow direction and the tangent to the 
parabolic surface and a is the angle, lying on the x-z plane, 
between the z-axis and the line of sight. This latter can be 
inferred by measuring in the radial velocity map (Figure [T4l 
left panel) the projection D' over the sky plane of the dis- 
tance between the line emission peak and the maximum 
radial velocity (D' =10±2 pixels). We have: 



D' = Dcos(x -a) = \A4 Q2 . + 2z max b x cos(x - a) (A.l) 

where z max is the distance along the z-axis at which 
the radial velocity reaches its maximum. We firstly esti- 
mated the b parameter. Since the molecular hydrogen does 
not emit over the entire bow surface, there is a leading 
edge which divides the bow into two different zones: a 
dissociation cap, beyond which the emission comes from 
atomic/ionic elements and a molecular hydrogen emitting 
flank. Therefore, it is possible to define as Zdis the position 
along the z-axis where the leading edge plane is located: 
certainly it depends on the shock velocity at the bow apex, 
Vbow , on the H2 breakdown velocity, Vdis and on the shape 



of the bow surface. Setting Vdis = 
for a generic angle 8, tan0=dR/dz: 



Vdi. 



I'bo 



dR/dz 



v bow sin0 dis and, since 



(A.2) 



Substituting the parabolic equation in this formula, we 



get: 



'bow 



(A.3) 
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and 



b = 



dis 




(A.4) 
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In ea lA.41 Rdis is the radius of the cap beyond which H2 
dissociates, that we have measured both from the 2.122/mi 
(Figure [T5|) and 1.644/im emission contours fFigurc [T4")) . 
Both of them give Rdis ~ 20 pixels. We also have esti- 
mated Vbow = 115±5kms _1 from the FWZI of the 1.644/im 
(Hartigan, Raymond & Hartmann, 1987). Substituting 
these values in eci lA.4[ we b as a function of v<ns only, that 
we have taken as a free parameter. 

The above expression of b is then substituted in the 
parabolic equation to derive z max : 



Following Hartigan, Raymond & Hartmann (1987) and 
considering that the bow is seen from the back (D99), we 
get Qmax = tt/2 - a/2. Thus: 



In conclusion, from equations lA.4l and lA.6l and by mea- 
suring x from simple geometrical considerations, we are 
able to express in eq lA.ll D' as a function only of a and 
Vdisi that have been routinely varied to match our mea- 
surement of D', as shown in Figure fl~5l 




(A.5) 




(A.6) 



